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Abstract 

We have recently proposed a Bethe Ansatz solution of the open spin-1/2 XXZ quan- 
tum spin chain with general integrable boundary terms (containing six free boundary 
parameters) at roots of unity. We use this solution, together with an appropriate string 
hypothesis, to compute the boundary energy of the chain in the thermodynamic limit. 



1 Introduction 



There has been considerable interest in the open spin-1/2 XXZ quantum spin chain with 
general integrable boundary terms lU |2] , whose Hamiltonian can be written as 



+ - sinh r] ^ coth a_ tanh /?_cr^ + cosech a_ sech ( cosh O^a^ + i sinh ] 



'1.1] 



— coth a+ tanh l3+a% + cosech a+ sech j3+{ cosh 6+a% + i sinh ^^+cr^) 
where t] is the bulk anisotropy parameter, and a± , f3± , 6± are free boundary parameters. § 

Except for the special case a± oi (3± oo when the boundary terms become diagonal 
[31111 [5], the boundary terms break the bulk U{1) symmetry generated by S*^; i.e., the model 
has no continuous symmetry. For generic values of boundary parameters, this model does 
not seem to have a simple pseudovacuum, which precludes constructing a conventional alge- 
braic Bethe Ansatz solution. Being associated with the spin-1/2 representation of f/g(sM(2)), 
this model is but the simplest of an infinite hierarchy of more complicated integrable quan- 
tum spin chains involving higher- dimensional representations and/or higher-rank algebras. 
Hence, solving the former model is presumably a prerequisite for solving any of the latter 
ones. This model also has numerous applications in statistical mechanics, condensed matter 
and quantum field theory. 

A Bethe Ansatz solution of this model was found in [6]- [9] for the case that the boundary 
parameters obey the constraint 

1-62 

a_ + ti(3_ + t2a+ + e'i(3+ = to{6_ - 9+) + r]k ^ ^ — m mod (2i7r) , eie2e3 = +1 , (1.2) 

where = ±1, and k is an integer such that \k\ < N — 1 and A^ — 1 + fc is even. Completeness 
of this solution is not straightforward, as two sets of Bethe Ansatz equations are generally 
needed in order to obtain all 2^ levels [8]. Related work includes [T0]-|16]. 

There remained the problem of solving the model (11. ip when the constraint (II. 2p is not 
satisfied. Building on earlier work [171 [IB], we recently proposed in [19] a solution of the 
model for arbitrary values of the boundary parameters, provided that the bulk anisotropy 
parameter has values 



^Under a global spin rotation about the z axis, the bulk terms remain invariant, and the boundary 
parameters 9± become shifted by the same constant, 0^ ^ 9± -\- const. Hence, the energy (and in fact, the 
transfer matrix eigenvalues) depend on 6± only through the difference 9- — 9+. 
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where p is a positive integer. Hence, g = is a root of unity, satisfying q^^^ = — 1. 

As is well known, for both the closed chain and the open chain with diagonal boundary 
terms, the eigenvalues of the Hamiltonian (and more generally, the transfer matrix) can 
be expressed in terms of zeros ("Bethe roots") of a single function Q{u). This is in sharp 
contrast with the solution [19], which involves multiple Q functions, and therefore, multiple 
sets of Bethe roots. The number of such Q functions depends on the value of p. (Generalized 
T — Q equations involving two such Q functions first arose in [18] for special values of the 
boundary parameters.) 

The solution [19] has additional properties which distinguish it from typical Bethe Ansatz 
solutions: the Q functions also have normalization constants which must be determined; and 
the Bethe Ansatz equations have a nonconventional form. Given the unusual nature of this 
solution, one can justifiably wonder whether it provides a practical means of computing 
properties of the chain in the thermodynamic {N oo) limit. To address this question, we 
set out to compute the so-called boundary or surface energy (i.e., the order 1 contribution to 
the ground-state energy), which is perhaps the most accessible boundary-dependent quantity. 
For the case of diagonal boundary terms, this quantity was first computed numerically in 
[1], and then analytically in [20] . 

We find that the boundary energy computation is indeed feasible. The key point is 
that, when the boundary parameters are in some suitable domain, the ground-state Bethe 
roots appear to follow certain remarkable patterns. By assuming the strict validity of these 
patterns ("string hypothesis"), the Bethe equations reduce to a conventional form. Hence, 
standard techniques can then be used to complete the computation. We find that our final 
result (13.281) for the boundary energy coincides with the result obtained in [21] for the case 
that the boundary parameters obey the constraint (II. 2p . and in [22] for special values [T7l[T8] 
of the boundary parameters at roots of unity. 

The outline of this paper is as follows. In Section [2|, we briefly review the Bethe Ansatz 
solution [19] of the model ( 11. ip at roots of unity ( 11. 3p . In Section [3l we treat the case of even 
p, followed by the case of odd p in Section HI This is followed by discussion of our results 
and a brief outline of our future work in Section [5l 

2 Bethe Ansatz 

In this section, we briefly recall the Bethe Ansatz solution |T9|. In order to ensure hermiticity 
of the Hamiltonian (1.1), we take the boundary parameters (3± real; a± imaginary; 9± 
imaginary. We begin by introducing the Ansatz for the various Q{u) functions that appear 
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in our solution, which we denote as aj{u) and bj{u): 

aj{u) = Aj sinh(M — u^j^'^) , bj{u) = Bj JJ^ sinh(u — u^^^^) , 

k=l k=l 

J = l,...,[^^\, (2.1) 

where {u^^^^ , uf'^} are the zeros of aj{u) and bj{u) respectively, and [ J denotes integer part. 
If p is even, then there is one additional set of functions corresponding to j = | + 1, 

(ap ) (ap ) 

aE+i(n) = v4£+i JJ sinh(M - ) sinh(u + ^"^^ ) , 

k=l 



&£+i(m) = Be+i Yl sinh(M - ) sinh(M + ) . (2.2) 

k=l 

The normalization constants {y4j,i?j} are yet to be determinedly. We assume that is 
even, in which case the integers Ma , Mf, are given by 

Ma = - + 2p, M, = -+p-l, (2.3) 

It is clear from (12.11) . (12. 2p that aj{u) and bj{u) have the following periodicity and crossing 
properties, 

p 

aj{u + in) = aj{u), bj{u + in) = bj{u) , j = 1 , . . . , [-J + 1 , (2.4) 

aE+i(-M) = aE+i{u), bE+i{~u) = bE+i{u) . (2.5) 

The zeros of the functions {aj{u)} and {bj{u)} satisfy the following Bethe Ansatz equa- 
tions 



h,{-ur'-ri) ^ a,{-ur) + gijun Yjurr b^j-un 



(2.6) 



j = 2,... ,L^J+1, (2.7) 



and 



+ (J - l)r/) a.MrV 2- 



(2.8) 



^^"^''-^"^ - ^-^("^') , = 2,...,L^J+1, (2.9) 



h{uf^H{j-l)v) &,+i(«^V '■■■ 2- 



^One of these normalization constants can be set to unity. 
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where ap^2{'^) = ap{~u) and a^+s (m) = ap+i{—u) for even and odd values of p, respectively, 

2 2 2 2 

and similarly for the 6's. Moreover, 

h{u) = hoiu) hi{u) , (2.10) 
where Hq^u) and hi{u) are as follows 

ho{u) = smh^^{u + ri)- 



27V/ smh.{2u + 21]) 



sinh(2M + rj) ' 

hi{u) = —4 sinh(u + a;_) cosh(u + sinh(M + a+) cosh(M + /?+) . (2.11) 
We also define the quantities 

^i(n) = 2sinh(2(p + l)u) (2.12) 

and 

2 

Y{uf = J2f^k cosh'' {2{p+l)u) . (2.13) 

fc=0 

Explicit expressions for the coefficients fik in (12.131) . which depend on the boundary param- 
eters, as well as for the function /i(m), are listed in the Appendix for both even and odd 
values of p. 

Moreover, there are additional Bethe-Ansatz-like equations 

= «2(-|), (2.14) 

a,-iiil~j)v) = «,+i((^-J», j = 2,...,l|J+1, (2.15) 
which relate the normalization constants {^j}; and also 

= &2(-|), (2.16) 

bj-i{{l-j)v) = J = 2,...,l|J+1, (2.17) 

which relate the normalization constants {Bj}. There are also equations that relate the 
normalization constants Ai and Bi, such as 

/i(-a_ - T]) bi{a- +r]) = -gi{-a- - r]) ai(a_ + r]) . (2.18) 
The energy eigenvalues of the Hamiltonian (11.11) are given by 

E = - sinh 7] ^ [coth(Mf + (j - l)r]) - coth(Mf ^"'^ + (j - 1)7]) 



J = 2,...,L^J, (2.19) 



where Eq is defined as 

Eq = - sinh?7 (cotha_ + tanh + cotha+ + tanh /?+) + -{N — 1) cosh 77 . (2.20) 
For even p, there is one more expression for the energy corresponding to j = | + 1, 

E = -smhrji coth(M; H ) — cotli^w^ ^ ) 

2 1 2 2 

- X]coth(«;^ +y)|+^o. (2.21) 

There are also similar expressions for the energy in terms of a roots {u^""^^} [19j. 

3 Even p 

In this section, we consider the case where the bulk anisotropy parameter assumes the values 
fll.Sp with p even, i.e., rj = y , y , . . .. We have studied the Bethe roots corresponding to 
the ground state numerically for small values of p and along the lines of |8]. We have 
found that, when the boundary parameters are in some suitable domain (which we discuss 
further below Eq. (13.281) ). the ground state Bethe roots {u^k^\uf^^} have a remarkable 
pattern. An example with p = 2 , = 4 is shown in Figure [H Specifically, these roots can 
be categorized into "sea" roots, {v'^^'^'\v^^''^^} (the number of which depends on A^) and the 
remaining "extra" roots, {wf^"'^'''\w'^^''^^} (the number of which depends on p) according to 
the following pattern which we adopt as our "string hypothesis" . 

3.1 Sea roots {vt^'''\v^^^'^} 

Sea roots of all {aj{u) , bj{u)} functions for any even p are summarized below, 

j = l,...,| + l, (3.1) 

where Vk are real and positive. In Figure [H the sea roots are indicated with red stars. 

Note that the real parts (i-u^) are independent of j. This, as we shall see, greatly sim- 
plifies the analysis. Furthermore, for each sea root with real part +Vk, there is an additional 
"mirror" sea root with real part —Vk, for a total of A^ sea roots, provided j 7^ §"'"-'-■ 
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Figure 1: Ground-state Bethe roots ioi p = 2, N = 4, «_ = 0.604i, a+ = 0.535i, 
/?_ = -1.882, p+ = 1.878, e_ = 0.6i, 9+ = 0.7i. 

j = I + 1, there are only y sea roots +Vk + y (i.e., just the root with positive real part) 
due to the crossing symmetry (12. 5p of the functions a£+i(u) and bp^i{u). % 



3.2 Extra roots {^^^"""'^ , w^^^'^} 

We next describe the remaining extra Bethe roots for even p, the number of which depends 
on the value of p. In Figure [H the extra roots are indicated with black circles. Since the 
functions aj{u) and bj{u) have a different number of such extra roots, we present them 
separately. The extra roots of the bj{u) functions have the form 



w 



2p 



2k 



V: 



k = 1 



,p-l 



p 



(3,2) 



■^Hence, strictly speaking, we should write the j = f + 1 equation in (|3.ip separately, keeping only the + 
roots. However, in order to avoid doubling the number of equations, we commit this abuse of notation here 
and throughout this section. 
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The real parts of the roots, Wk, are not all independent. Instead, they are related to each 
other pairwise as follows, 

P 



Wk = Wp-k , 



k = l - - 1 

2 



(3.3) 



Only We remains unpaired. This property proves to be crucial for the boundary energy 
calculation. 

There are two types of extra roots of the aj{u) functions: 



w 



±(a,,l) 



Wk ^"^k + 



2p+l-2k 



r], k = 1,... ,p-l 



w 



±(aj,2) 



±Wo + 



2p + 3-2k 



P 



7], k = 1,... 



(3.4) 



Note that the extra roots of the first type {w'^^"'^'^^} coincide with the b roots {w^^^^^}; and 
that the extra roots of the second type {w^''"^'^'*} form a l)-string", with real part Wq. 

As previously remarked, for j = | + 1, only the roots with the + sign appear. 



3.3 Boundary energy 



We now proceed to compute the boundary energy. Using the expression (12.191) for the energy 
and our string hypothesis, we obtain (for p > 2) 

JV 

1 r ^ r 

E = - sinhril ^ ^coth(t;J^*^^ + (j - l)r/) + coth(t;;^*^^ + (j - l)r/) 

k=l 

- coth(t;+^''^-^^ + (j - l)r]) - coth(t;;(''^-^^ + (j - l)r]) 
+ J2 [coth(«;+^''^^ + (j - l)r]) + coth{w~^'''^ + (j - l)r/) 



k=l 



- coth(u;+^''^-'^ + (j - 1)t]) - coth(u7^ + (j - l)r]) 



(bj-i) 



3 = 2,... ^. 



Recalling (13. ip and (13. 2p . this expression for the energy reduces to 

1 



JV 
2 



E = sinh^ 11 y — 



^ sinh(t;fc - 2) sinh(t;fc + 2) 



+ £"0 , Vk>^. 



(3.5) 



(3.6) 
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independently of the value of j. Since the extra roots wf^'' are independent of j, their 
contribution to the energy evidently cancels, leaving only the sea-root terms in (13.61) . The 
same result can also be obtained (for p > 2) from the energy expression (I2.2ip . 

We turn now to the Bethe Ansatz equations, on which we must also impose our string 
hypothesis. Choosing j = | + 1 in (12.91) with u] equal to the sea root ^ = f / + y, 
we obtain 



where we have made use of the fact 6£+2('^^) = bt^—u). More explicitly, this equation reads 
'smh{vi + sinh{2vi + r]) smh{vi — f + a;„) cosh({;; — ^ + 



sinh(t>/ — f )/ sinh(2-y/ — rj) sinh('y/ + ^ — «-) cosh(fi + ^ — /?-) 

JV 

sinh({;/ - f + a+) cosh('ii; - ^ + (3+) _ -q sinh(-yi - Vk + r]) sinh(ii; + Vk + r]) 



X 



sinh(t;/ + ^ - a+) cosh ({5; + 2 _ ^_^) 11 sinh(t;; - Vk - f]) sinh(fi + Vfc - ^) ' 

Z = l,---,y, ?}fc>0. (3.8) 

In obtaining this result, we have made use of the fact that the normalization constant Bv of 
the function &£(«) cancels, and also that the contribution from the extra roots on the RHS 
cancel as a consequence of the relation (13.31) among their real parts. 

Remarkably, as a consequence of our string hypothesis, our non- conventional Bethe 
Ansatz equations have reduced to a conventional system (13. 8p . which can be analyzed by 
standard methods. However, before proceeding further with this computation, it is worth 
noting that the same equations can also be obtained starting from any j > 1. To see this, 
we first observe that the {Aj} normalization constants are all equal, and similarly for the 
{-Bj} normalization constants, 

A^=A2 = ... = Az+i , Bi = B2 = ... = Be+, . (3.9) 



This result follows from the Bethe- Ansatz-like equations (12.14p - (l2.17p and the string hypoth- 
esis. For example, using (13. ip and (13. 2p in (12.160 . and remembering the relation (13. 3p among 
the real parts of the extra roots, we obtain Bi = B2. Hence, choosing in (12. 9p to be a 
sea root f^'^*'^^'' for any j G {2 , . . . , | + 1}, we again arrive at (13. 8p . Moreover, in view of the 
identity 

,-2 ^ + 1 (310) 
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where v^^"'^'^ = v^^^^^ is a sea root, the same result fl3.8l) can also be obtained from fl2.7l) . 

In the thermodynamic (A^ oo) limit, the number of sea roots becomes infinite. The 
distribution of the real parts of these roots {vk} can be represented by a density function, 
which is computed from the counting function. To this end, following [2T| [22] and references 
therein, we define some basic quantities 

^' sinh/i(A-f) ^ 2/i^ cosh/i(A-f) ^ ' 

which allow us to rewrite the Bethe Ansatz equations (13.81) in a more compact form. 



JV 
2 



i\\2N+i X e2a_-i(A/) e2a+-i(Ai) TT /' \ \ ^ m i \ ^ 

ei(AO ^ ^i(A,) — — - = - I I ealAz - Afc) e2(A, + Afe) , 

^i+2i6_(A/) gi+2ib+[M) ^Ji 

/ = l,---,y, (3.12) 

where we have set vi = /lA/, t] = ,a± = ifia± , /5± = fib±. Note that the parameters fi, a±, 
b± are all real. 

Taking the logarithm of (13.121) . we obtain the desired ground state counting function 
h(A) = 7^|(2Ar+ l)gi(A) +ri(A) + g2a_-i(A) -ri+2i6_(A) + g2a+-i(A) -ri+2i6+(A) 

JV 

- X][g2(A-A,) + g2(A + Afc)]}, (3.13) 

k=l 

where gn(A) and r„(A) are odd functions defined by 

g„(A) = TT + 2 lne„(A) = 2 tan~^ (cot(n/i/2) tanh(/xA)) , 

r„(A) = ^lnf7n(A). (3.14) 

Defining A_fc = — A^, we have 

JV JV 
2 2 

- E fe(A - Afc) + g2(A + Afc)] = - J2 92(A - A,) + g2(A) . (3.15) 



k=l 



The root density p(A) for the ground state is therefore given by 

1 r/h 1 r 

P(A) = -^ = 2a,{X)-j dX'a,{X-X')p{X') + -[a,{X) + W{X) 

+ a2(A) + a2a_-i(A) - &i+2i6_(A) + a2a+-i(A) - 6i+2i6+(A) , (3.16) 



*Only the first set of Bethe equations (|2.6p . (|2.8p do not seem to reduce to ([3 



9 



where we have ignored corrections of higher order in when passing from a sum to an 
integral, and we have introduced the notations T 



an. (A) 

&n(A) 



Id n sin(n/i) 

27r (iA " 77 cosh(2|uA) — cos(n/i) 



1 d 



r„(A) 



sin(?T,yu) 



(3.17) 



271 dX ^ TT cosh(2yuA) + cos(nyu) 
The solution of the linear integral equation (13.161) for p(A) is obtained by Fourier transforms 
and is given by T 



p(A) = 2s(A) + 1/?(A) 



where 



and 



1 f 



^11 

dcj e"*"^^ , , , , = — — 

2cosh(cj/2) 2cosh(7rA) 



(3.18) 
(3.19) 



(1 + 02(0;)) 

+ a2a_-l(t^) + a2a+-l(t^)} 



ai(c^) + 61 (c^) + 02(0;) - 6l+2*_(t^) - &l+2ib+(t^) 



with 



dn{uj) 



sgn(?T, 



sinh ((z/ — \n\)u!/2) 



bn{uj) 

where u = - = p + 1. 



sinh {i'uj/2) 
sinh {nuj/2) 



sinh {vuj/2) 



< < 2z/, 
< 3fJe n < . 



(3.20) 

(3.21) 
(3.22) 



Expressing the energy expression (13. 6p in terms of the newly defined quantities and letting 
become large, we obtain 



E 



27C sin fi 



N_ 

2 



TTSm/i 



^ai(Afc) + Eo 



k=l 



Tc sm /i 



JV 

2 



{ ^ ai(Afe) -ai(0)} +Eo 



fx 



N 



I dA ai(A) p(A) -ai(0)| + -(A^- l)cos/i 



+ - sin p (cot /ia_ + i tanh fib- + cot /ia+ + i tanh /i&+) , 



(3.23) 



''These new functions an(A) and bn{X) should not be confused with the Q functions aj{u) and bjiu) 
appearing earher. We apologize for this unfortunate coincidence of notations. 
^Our conventions are 



/(^) 



e'"^ /(A) d\.. 



ZTT 



^ — iujX 



/(w) dw . 
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where again we ignore corrections that are higher order in Substituting the result 

f l3.18p for the root density, we obtain 



E — Ebulk + Eboundary , (3.24) 

where the bulk (order A^) energy is given by 

2A^7rsin/i /'°° „ ,,,,,, 1.. 
Ebulk = / rfA ai(A) s(A) + -A/ cos/i 

= -iVsinV / dX — — + -iVcos/i, (3.25) 

j_oo [cosh(2/iA) — COS /ij cosh(7rA) 2 

which agrees with the well-known result [23]. The boundary (order 1) energy is given by 

Eboundary = " — COS /i H — siu /i (cot /ia_ + i tauh + cot /ia+ + i tanh fib^){^.26) 

where / is the integral 

d\ ai(A) [R{X) - 5(A)] = — / dcu ai(cu) R{lu) - 1 
duj s{uj)< di{uj) + bi(uj) — 1 



1 

2^ 



oo 

oo 



-bl+2ib_iuj) - 6l+2i5+(c^) + a2a_-i(t^) + a2a+-i(t^)| • (3.27) 

We further write the boundary energy as the sum of contributions from the left and right 
boundaries, Eboundary = -Sfo^^tndary + ^boundary ^hc energy contribution from each boundary 
is given by 



, ^ sin/i 1 f sinh((^-2)a;/4 ) 1 

boundary J 2 COsh(c^/2) I 2 sinh(z/C^/4) 



+ osn(2a^ ^^ sinh((z/- |2a±-l|)a;/2) ^ sinh((2z&± + 1)^/2) | 

sinh(i/ci;/2) sinh(i/u;/2) i 

+ - sin yU (cot yua± + z tanh /i6-|-) — - cos /i . (3.28) 



This result can be shown to coincide with previous results in 

We emphasize that the result (13.281) has been derived under the assumption that the 
Bethe roots for the ground state obey the string hypothesis, which is true only for suitable 
values of the boundary parameters. For example, the shaded areas in Figures [2] and [3] denote 
regions of parameter space for which the ground-state Bethe roots have the form described 
in Sections |3 . 1 1 and [3T2l The a± and [3± parameters are varied in the two figures, respectively. 
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Figure 2: Shaded area denotes region of the (Sm a;+ , S^m a_) plane for which the 
ground-state Bethe roots obey the string hypothesis for p = 2, = 2, 



/3- 



-1.882, /5+ = 1.878, = O.Qi, 9+ = 0.7i. 




3 
2.5 
2 

l.S 

1 

0.5 



Figure 3: Shaded area denotes region of the ,/?_) plane for which the ground- 
state Bethe roots obey the string hypothesis for p = 2, = 2, a_ = 
-1.818i, a+ = 2.959i, 9. = 0.7i, 9+ = 0.6i. 

4 Odd p 

In this section, we consider the case where the bulk anisotropy parameter assumes the 
values (11.31) with p odd, i.e., 77 = ^,^,.... As for the even p case, for suitable values of the 



boundary parameters, the ground state Bethe roots {ul 



} have a regular pattern. An 



example with p = 3 , A^ = 4 is shown in Figure HI As before, these roots can be categorized 
into sea roots (the number of which depends on A^) and extra roots (the number of which 
depends on p) according to the following pattern which we adopt as our "string hypothesis" . 



4.1 Sea roots {i;^^"^-*, t'^^^^'*} 



Sea roots of all {aj{u) , bj{u)} functions for odd p are given by 



N 
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Figure 4: Ground-state Bethe roots for p = 3, N = 4, a_ = 1.554i, a. 



0.948i, 



/3_ = -0.214, /3+ = 0.186, ^_ = 0.6i, ^^ 



0.72. 



p+ 1 



(4.1) 



where Vk are real and positive. In Figure HJ the sea roots are indicated with red stars. 

As in the even p case, the real parts (±?)fc) are independent of j. This again provides 
simplification to the analysis. In contrast to the even p case, now none of the functions 
{aj{u) , bj{u)} has crossing symmetry. Hence, there are N sea roots for all values of j. 



4.2 Extra roots {w^^^'''\wf^^} 

We now describe the extra Bethe roots for odd p. In Figure HJ the extra roots are indicated 
with black circles. We start with the p — 1 extra roots of the bj (u) functions: 



w 



w. 



p-i 



±Wk + {p-k)r], /c = l,...,p-2. 

> + 2-2j\ , p+1 



p-i 



rj. 



3 



(4.2) 
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Similarly to the even p case, the real parts of the extra roots are related to each other 
pairwise, 

Wk = Wp-fc-i , A; = 1 , . . . , , (4.3) 
so that only remains unpaired. 

2 

Similarly, the extra roots of the aj{u) functions are as follows, 
^±(aj,i) ^ w^''^^^ = ±Wk + {p - k) 1] , k = l,...,p-2, 

±(a„l) ±{bj) ,~ f p+2-2j\ 

Wp-i = = ±^P-i + ( 2 J ^' 

^±(a„2) ^ + fc = l,...,p+l, j = l^...^P±l. (4.4) 



As in the even p case, the extra roots of the first type {w^^""^'^^} coincide with the b roots 
{wf^'^^^}. Moreover, the extra roots of the second type {w^^"'^''^^} form a "(p + l)-string", 
with real part wq. 

However, in contrast to the even p case, some of the extra roots (namely, w^j!.i^ and 
w^p-i) depend on the value of j. Hence, as we shall see, these extra roots will not cancel from 
either the energy expression or the Bethe equations. Nevertheless, the contribution of these 
roots to the boundary energy will ultimately cancel. 

4.3 Boundary energy 

As in the case of even p, we use the energy expression 02.191) and the string hypothesis to 
obtain (for p > 3) 



1 r r 

E = - smhr]l J2 [coth(4^*'^ + (j - 1)^) + coth(w;^*^^ + (j - l)r]) 

^ k=l 

- coth(4^''^-^^ + (j - l)r]) - coth(w;^'^-^^ + {j - l)r]) 
+ J2 [coth(w;^^^^^ + (j - l)r]) + coth(M;^ ^^^^ + (j - l)r/) 



coth(u;+^''^-'^ + (j - 1)t]) - coth(u;^ ^'''^'^ + (j - l)r]) 



J = 2.....f^. (4.5) 
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Recalling (14. ip and (14. 2p . this expression for the energy reduces, independently of the value 
of j, to 

2 1 2 sinh^ r] 

E = smh r] > . , , „, . , ^ ; r— — r + -Eo , (4.6) 

^ sinh(ffc — ^) sinh(ffc + ^) cosh?] + cosh(2wp_i) 

where Vk ,Wp-i > 0. As already anticipated, the expression for the energy depends on the 
extra root Wp^i as well as on the sea roots. 

Turning now to the Bethe Ansatz equations, following similar arguments as for the even 
p case, we find again that the A normalization constants are all equal, and similarly for the 

5's, 

A^ = A2 = ... = Ap+i , Bi = B2 = ... = Bv+i . (4.7) 

2 2 

Choosing uf^^ in (12.90 to be a sea root f;^^^^^ for any j G {2 , . . . , ^^}, we obtain 

smh.{vi + ^)\^^ smh{2vi + r/) smh{vi - f + «-) cos\i{vi - ^ + /5-) 
sinh({)/ — f)/ sinh(2{)/ — rf) sinh({)/ + f ~ '^-) cosh({;; + f ~ /^-) 

sinh({;/ - f + a+) cosh.{vi - ^ + /?+) sinh({;/ - Wp^i - smh{vi + Wp^i - ^rf) 

sinh({;/ + ^ - a+) cosh.{vi + ^ - (3+) sm\i{vi - Wp-i + smh.{vi + Wp-i + ^r]) 

JV 

nsmh{vi - Vk + v) s^'ohivi + Vk + T]) N ^ ^ 

—rp z r—TT^—^ r, / = 1,---, — , Vk,Wp_i>0. (4.8) 
smh(t;, -Vk-rj) sm\i[vi + Vk - rj) 2 

In a compact form, this result can be written as 

ei(Ai) ^ ^i(Ai) — — - = - [ep_i(A; - A) ep_i(Ai + A)J 

JV 

ne2(A/-Afc) e2(A, + Afc), / = 1, (4.9) 



X 

fc=i 



2 

where Wp_i = /iA. The corresponding ground state counting function is given by 

h(A) = ^|(2A^ + l)?i(A) + ri(A) + g2a_-i(A) - n+2ib-{\) + g2a+-i(A) - n+2^b+{>^) 



N 
2 



+ gp_i(A - A) + gp_i(A + A) - ^ [g2(A - A,) + g2(A + A,)] } . (4.10) 

k=l 

Following similar procedure as before, we arrive at the root density for the ground state 
p(A) = 2ai(A) - r d\' a2(A - A') p(A') + ^ L (A) + 6i(A) + a2(A) (4.11) 

+ a2a--i(A) - 6i+2ib_(A) + a2a+-i(A) - &i+2ib+(A) + ap-i(A - A) + ap_i(A + A) 
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where as before higher order corrections in are ignored when passing from a sum to an 
integral. This yields 



where now 



R{uj) 



1 



p(A) = 2s(A) + ^i?(A) , 



(4.12) 



(1 + 02(0;)) 

+ a2a_-i(t^) + a2a+-i(c^) + 2 cos(Au;) ap-i(u;) [> . 



(4.13) 



The energy expression ( 14. 6p yields, as iV — > oo, 

iV 

2'^^'"^.{f:.,(A.)+6,(A)} + £„ 



TT Sm yU 



fc=l 



2 



J] ai(Afe)-ai(0) + 26i(A)}+Eo 



= -!^^|Ar r dX ai(A) p(A) - ai(0) + 26i(A)| + ^(A^ - 1) cos/i 

+ - sin /i (cot /ia_ + i tanh yu6_ + cot /ia+ + i tanh /i6+) . 
2 



Substituting (14.121) for the root density, we again obtain 



(4.14) 



(4.15) 



where the bulk (order A^) energy is again given by fl3.25p . The boundary energy is now given 
by 



E 



boundary 



TT sin 1 1 . • 1 , • 1 , N 

1 COS /i H — sm p (cot /ia_ + t tanh /iO_ + cot /ia+ + t tanh /iO+) 

1^ 2 2 



27r sin /i 



&i(A), 



(4.16) 



where I is now the integral 

/ = / c/A ai(A) [i?(A) -5(A)] 



1 /"°° 

— / cico' ai(u;) R{uj) — 1 
-bi+2ib.{uj) - bi+2ib+{uj) + a2a_-i{u}) + a2a+-i{uj) + 2cos(Alj) dp_i{uj) \ . (4.17) 



1 

2^ 
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Using the fact that s{uj)dp_i{uj) = —bi{uj), we see that there is a perfect cancellation of the 
last term in fl4.16p which depends on the extra root A. Thus, as in the even p case, there 
is no contribution to the boundary energy from extra roots. Proceeding as before, we find 
that the energy contribution from each boundary is again given by fl3.28p . thus coinciding 
with previous results in ^ll 122] . 

As for even p, the derivation here is based on the string hypothesis for the ground-state 
Bethe roots, which is true only for suitable values of boundary parameters. For example, 
the shaded areas in Figures and denote the regions of parameter space for which the 
ground-state Bethe roots have the form described in Sections 14.11 and 14. 2[ The a± and /3± 
parameters are varied in the two figures, respectively. 





1 2 



Figure 5: Shaded area denotes region of the {Qma+ ,Qma-) plane for which the 
ground-state Bethe roots obey the string hypothesis for p = 3, N = 2, 
p_ = -0.85, P+ = 0.9, 9_ = 0.6'i, 9+ = 0.7i. 




Figure 6: Shaded area denotes region of the ,/?_) plane for which the ground- 
state Bethe roots obey the string hypothesis for p = 3, = 2, a_ = 1.2i, 
a+ = 0.98i, = 0.7i, 9+ = 0.6i. 



17 



5 Discussion 



We have studied the ground state of the general integrable open XXZ spin-1/2 chain f ll.ip in 
the thermodynamic hmit, utihzing the solution we found recently in jl9j. In contrast to the 
earlier solution [6]- [9], this solution does not assume any restrictions or constraints among 
the boundary parameters. However, the bulk parameter is restricted to values corresponding 
to roots of unity (11.31) . The key to working with this solution is formulating an appropriate 
string hypothesis, which leads to a reduction of the Bethe Ansatz equations to a conventional 
form. While the idea of using a string hypothesis to simplify the analysis of Bethe equations 
is as old as the Bethe Ansatz itself, the particular patterns appearing here are perhaps 
unparalleled in their rich structure. 

The boundary energy result f l3.28p was obtained previously [2T] for bulk and boundary 
parameters that are unconstrained and constrained, respectively; and we have now obtained 
the same result for the reversed situation, namely, for bulk and boundary parameters that 
are constrained and unconstrained, respectively. Hence, this result presumably holds when 
both the bulk and boundary parameters are unconstrained (within some suitable domains). 
Indeed, for the boundary sine-Gordon model [2], which is closely related to the open XXZ 
chain, the expression [23] for the boundary energy is valid for general values of the bulk and 
boundary parameters. In view of the spectral equivalence between systems with diagonal 
and nondiagonal boundary interactions noted in [HI [ISl [16], it may be interesting to try to 
relate our boundary energy result with the corresponding result [20] for diagonal boundary 
interactions. 

Having demonstrated the practicality of this solution, we now expect that it should be 
possible to use a similar approach to analyze further properties of the model, such as the 
Casimir energy (order 1/A^ correction to the ground state energy), and bulk and boundary 
excited states. 

There is an evident redundancy in the solution which we have used here: there are many 
equivalent expressions for the energy (see, e.g., (12.191) . (12.211) ). and we find that the Bethe 
Ansatz equations (12.71) . (12.91) all become equivalent upon imposing the string hypothesis. 
Moreover, while there are various "extra" Bethe roots describing the ground state, they 
ultimately do not contribute to the boundary energy. All of this suggests that it may be 
possible to find a simpler and more economical solution of the model involving fewer Q 
functions. Ideally, one would like to find a solution for which neither bulk nor boundary 
parameters are constrained. 
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A Appendix 

We hst here exphcit expressions for the function fi{u) and the coefficients fik appearing in 
the text. We remind the reader that we assume throughout that is even. 

For even values of p, the function /i (m) appearing in the Bethe Ansatz equations (12. 6p - 
dZSD, (|2A8|) is given by 

Mu) = -2^-2^ ( 

sinh {{p + l)a_) cosh {{p + l)/5_) sinh {{p + l)a+) cosh {{p + l)/3+) cosh^ {{p + l)u) 

— cosh {{p + l)a_) sinh {{p + l)/3-) cosh {{p + l)a+) sinh {{p + l)/3^) sinh^ {{p + l)u) 

— cosh((p + l)(e„-e+))sinh^((p+l)M)cosh^((p + l)u)). (A.l) 

For odd values of p, 
h{u) = -2'-'^[ 

cosh ((p + l)a_) cosh ({p + l)/3_) cosh ((p + l)a+) cosh {(p + l)P+) sinh^ ((p + l)u) 

— sinh ((p + 1)q;-) sinh {{p + l)/3_) sinh {{p + l)a+) sinh {{p + l)/5+) cosh^ {{p + l)u) 
+ cosh((p+l)(0_-^+))sinh2((p+l)M)cosh2((p+l)u)). (A.2) 

For both even and odd values of p, these functions have the properties 

f,{u + v) = fiiu), f,(-u) = Mu). (A.3) 

The coefficients fik appearing in the function Y{u) fl2.13p are given as follows for even 
(upper sign) and odd (lower sign) values of p. 

fio = 2-'pS^-l-cosh\{p + l){9_-9+)) 

- cosh(2(p + l)a_) cosh(2(p + l)a+) =f cosh(2(p + 1)«_) cosh(2(]9 + l)/?_) 
=F cosh(2(p + l)a+) cosh(2(p + l)/3_) =f cosh(2(p + l)a_) cosh(2(p + l)/3+) 
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T cosh(2(p + 1)«+) cosh(2(p + - cosh(2(p + !)/?_) cosh(2(p + 1)/?+) 



+ 



cosh((p + + a+)) cosh((p + - 



± cosh((p + - a+)) cosh((p + l)(/3_ + 



T 2cosh((p+l)(0_ -0+)) cosh((p + l)(a_ -a+))cosh((p+l)(^_ 



± cosh((p + l)(a_ + a+)) cosh((p + !)(/?_ + 



/ii = 2^-^^'|cosh((p+ l)(a_ cosh((p + ± a+)) 



+ cosh((p + ± cosh((p + 1){9_ - ^+)) 

T cosh((p + T /?+)) [cosh((p + ± 

+ cosh((p + l)(a_ ± a+)) cosh((p + 1)(^_ - 9+)) 

/i2 = 2-''Psmh\{p + l){e^-e+)). 



(A.4) 
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